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Abstract 

We study several azimuthal asymmetries in semi-inclusive deep inelastic scattering 
and in Drell-Yan, interpreting them within the formalism of the quark correlator, 
with a particular reference to T-odd functions. The correlator contains an 
undetermined energy scale, which we fix on the basis of a simple and rather general 
argument. We find a different value than the one assumed in previous treatments of 
T-odd functions. This implies different predictions on the Q 2 dependence of the 
above mentioned asymmetries. Our theoretical result on unpolarized Drell-Yan is 
compared with available data. Predictions on other azimuthal asymmetries could be 
tested against yields of planned experiments of Drell-Yan and semi-inclusive deep 

inelastic scattering. 



PACS numbers: 13.85. Qk, 13.88. +e 



1 Introduction 



Azimuthal asymmetries of remarkable size have been observed in various high energy- 
inclusive reactions, especially in unpolarized Drell-Yan[l, 2, 3], in singly polarized 
semi-inclusive deep inelastic scattering [4- 10] (SIDIS) and in inclusive production of 
(anti-)hyperons[ll-15] and pions[16-19] from singly polarized hadronic collisions. The 
interpretation of such asymmetries from basic principles of QCD is quite challeng- 
ing and has stimulated the interest of high energy physicists. In particular, in the 
present paper we focus our attention on the SIDIS and Drell-Yan asymmetries, which 
are somewhat analogous, since the two reactions are kinematically isomorphic. The 
theoretical activity about this subject is quite intense and lively, as witnessed by the 
numerous articles dedicated to the topic[20-70] in the last 15 years. 

An important element in the interpretation of such effects is the intrinsic trans- 
verse momentum of partons inside a hadron, whose crucial role in high energy reac- 
tions has been widely illustrated in the last years[23-25, 71-74]. Indeed, the transverse 
momentum is connected to the T-odd quark densities[20, 21, 27, 31-34], which pro- 
vide a quite natural interpretation of the above mentioned asymmetries[27, 31-34, 
42]. At the same time, the T-odd functions involve predictions of further azimuthal 
asymmetries in unpolarized and singly polarized inclusive reactions [29, 26]. 

These functions explain simultaneously [42] the remarkable cos2<f) asymmetry and 
the negligible coscj) Fourier component exhibited by unpolarized Drell-Yan data[l, 2, 
3], where is the usual azimuthal angle adopted in the phenomenological fits[l, 2, 3]. 
The term cos2(p may be just interpreted as a signature[42] of the pair of chiral-odd 
(and T-odd) functions involved in this picture. However, the current treatment of the 
T-odd functions does not reproduce the dependence of this asymmetry on the effective 
mass of the Drell-Yan lepton pair. More generally, some doubts have been cast on 
the Q 2 dependence of the transverse momentum distribution functions[75, 76, 77, 78], 
where Q is the QCD hard scale. This imposes a revision of the parameterization of 
the transverse momentum quark correlator, a fundamental theoretical tool for cross 
section calculations at high energies. This quantity - originally introduced by Ralston 
and Soper in 1979[79] and successively improved by Mulders and Tangerman[25, 23, 
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24] (see also the more recent contributions on the subject[30, 80, 81]) - consists of 
a 4 x 4 matrix. Therefore it may be parameterized according to the components 
of the Dirac algebra, taking into account the available vectors and hermiticity and 
Lorentz and parity invariance. The parameterization - whose coefficients are the 
quark distribution functions inside the hadron - includes an undetermined energy 
scale, Ho [73], usually assumed [23, 24, 25] equal to the mass of the hadron related to 
the active quark. We shall see that this choice is not unique, perhaps not the most 
appropriate in normalizing some "leading twist" functions. Alternatively, we propose 
A*o — Q/2, which explains quite naturally the Q 2 dependence of the unpolarized Drell- 
Yan asymmetry. Moreover, concerning the SIDIS and other Drell-Yan azimuthal 
asymmetries, we get predictions which contrast with those given by previous authors, 
and which could be tested against present[4, 5, 6, 7, 9], forthcoming [10] and future[82, 
83, 84, 85] data. 

Here we shall not study all azimuthal asymmetries considered in the literature [86, 
87, 88, 89], we shall limit ourselves to SIDIS of unpolarized or longitudinally polarized 
lepton beams off unpolarized or transversely polarized targets, and to unpolarized or 
singly polarized Drell-Yan, with transverse polarization; moreover, we shall consider 
just the asymmetries usually classified as leading twist [25, 27, 29]. 

The paper is organized as follows. In sect. 2 we give the general formulae for 
the SIDIS and Drell-Yan cross sections, introducing the formalism of the correlator; 
in particular we illustrate in detail the T-odd functions. Sect. 3 is dedicated to 
the theoretical formulae for azimuthal asymmetries. In sect. 4 we determine the 
parameter /i , by comparing the correlator with the quark density matrix in QCD 
parton model. Such a determination leads to predictions on the Q 2 dependence of 
the asymmetries, which we illustrate in sect. 5. In sect. 6 we compare our results 
with experimental data, as regards unpolarized Drell-Yan. Lastly we draw a short 
conclusion in sect. 7. 
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2 SIDIS and Drell-Yan Cross Sections 



2.1 General formulae 



Consider the SIDIS and the Drell-Yan reactions, i. e., 



lh A -> l'h B X 



and 



h A h B 



rrx, 



(i) 



where the Vs are charged leptons and the /i's are hadrons. Incidentally, these two 
reactions are topologically equivalent [22]. At not too high energies one can adopt 
one-photon exchange approximation, where the cross sections for such reactions have 
an expression of the type 

da _ (ATra) 2 

Here dT is the phase space element, a the fine structure constant and T = 



y (Pi ' P2) 2 — mlm 2 . the flux factor, p { and rrii (i — 1, 2) being the 4- momenta and 
the masses of the initial particles. Moreover L^ v and W^ v are respectively the leptonic 
and hadronic tensor. In particular, we have 



where I and J are the four-momenta of the initial and final lepton (in SIDIS) or 
of the two final leptons (in Drell-Yan). As regards the hadronic tensor, one often 
adopts, in the framework of the factorization theorem[90, 91, 22, 38], the so-called 
"handbag" approximation, where all information concerning the "soft" functions of 
the quark inside the hadrons is encoded in a parameterization of the quark-quark 
correlator, according to the various Dirac components [23, 25]. In this approximation 
the hadronic tensor reads 



Here c is due to color degree of freedom, c = 1 for SIDIS and c = 1/3 for Drell- 
Yan. §a and are correlators, relating the active (anti-)quarks to the (initial 
or final) hadrons h A and Kb- a, and b are the flavors of the active partons, with 
a = u,d, s,u, d,s and b = a in SIDIS, b = a in Drell-Yan; e a is the fractional charge 
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(3) 




a 




of flavor a. In Drell-Yan <&a and encode information on (anti-)quark distributions 
inside the initial hadrons: the x's are the longitudinal fractional momenta of the active 
quark and antiquark, is the transverse momentum of the active parton of h\ and 
q_i_ is the transverse momentum of the lepton pair. In SIDIS is replaced by the 
fragmentation correlator A[z, z(q± — p_i_)], describing the fragmentation of the struck 
quark into the final hadron h B (see subsect. 2.4). Here z is the longitudinal fractional 
momentum of hs with respect to the fragmenting quark and zq± is the transverse 
momentum of hg with respect to the virtual photon momentum. Approximation (4) 
holds for the hadronic tensor under the condition[42, 35, 36, 37] 

q± « Q, (5) 

where q± = |qj_|. Moreover we neglect the Sudakov suppression [92, 93], as usually 
assumed at moderate Q 2 [42, 54-62]. 



2.2 Parameterization of the Correlator 

The correlator for a nucleon may be parameterized according to the Dirac algebra, 
taking into account hermiticity and Lorentz and parity invariance. It is conveniently 
split into a T-even and a T-odd term, i. e., 

$ = $ e + $o, (6) 

where $ e is even under time reversal and $ D is odd under the same transformation. 
At leading twist one has[25, 94, 80, 81] 

V f I \ 

$ e ^ ^= + (A#il + A ± #i T )75/i+ + -h 1T "f5[p±, V>+\ I 

+ ^(^i L + \±hi T ) l5 [fi ± ,fi + l (7) 

$ ° " + (8) 

In formulae (7) and (8) we have used the notations of refs.[25, 23] for the "soft" 
functions*. n± are lightlike, dimensionless vectors, such that n + • n_ = 1 and whose 

*The correlator (6) has a different normalization than in rcf.[25], in accord with the definition 
(4) of the hadronic tensor. 
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space components are along (+) or opposite to (-) the nucleon momentum. Moreover 
the Pauli-Lubanski vector of the nucleon, denoted as S and such that S 2 = — 1, may 
be decomposed as 

S = + S ± . (9) 

Here P is the nucleon four-momentum, with P 2 = M 2 , A = —S-n and n = (0, 0, 0, 1) 
in the nucleon rest frame, taking the z-axis along the nucleon momentum. Thirdly, 

V = ^P-n-, \± = S- V±: (10) 

V± = P±/fJ>o, P± = P ~ ip ■ nJ)n + - (p ■ n+)n_ (11) 

and p is the quark four-momentum. Notice that the parameter V is similar to the one 
introduced by Jaffe and Ji[95, 96] with the same notation: denoting that parameter 
by Vjj, one has V = xVjj/V2. Lastly, the energy scale /i , encoded in the dimen- 
sionless vector rj±, has been introduced in such a way that all functions involved in 
the parameterization of $ have the dimensions of a probability density. This scale 
- defined for the first time in ref. [73], where it was denoted by mo - determines the 
normalization of the functions which depend on r]±; therefore /xq has to be chosen in 
such a way that these functions may be interpreted just as probability densities. We 
shall see in sect. 4 that taking /i equal to the rest mass of the hadron, as usually 
done[23, 24, 25], is not, perhaps, the most appropriate in this sense. Two observa- 
tions are in order about /i . First of all, it is washed out by integration over pj_ of the 
correlator, therefore it does not influence the common [95, 96] distribution functions. 
Secondly, we can reasonably assume that, for sufficiently large Q 2 , this parameter 
is independent of the perturbative interactions among partons: as we shall see, this 
conjecture can be proved rigorously. 



2.3 T-odd functions 

As explained in the introduction, the T-odd functions deserve especial attention. In 
particular, the two functions introduced in formula (8) may be interpreted as quark 
densities: corresponds to the quark transversity in an unpolarized (or spinless) 
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hadron, while is the density of unpolarized quarks inside a transversely polarized 
spinning hadron[27, 77]. 

A possible mechanism for generating these effects has been analyzed in detail, 
from different points of view, by various authors [35-41, 44-53]. In particular, the 
function f^, known as the Sivers function, may give rise to a single spin asymmetry, 
as predicted for the first time many years ago by Sivers[20, 21] as a consequence of 
coherence among partons. Essential ingredients for producing the effect are[35, 36, 37] 

a) two amplitudes with different quark helicities and different components (AL Z = 
1) of the orbital angular momentum; 

b) a phase difference between such amplitudes, caused, for example, by one gluon 
exchange between the spectator partons and the active quark, either before or after the 
hard scattering: owing to the different orbital angular momenta, the gluon interaction 
causes a different phase shift in the two amplitudes. 

Incidentally, a AL Z = 1 is connected to the anomalous magnetic moment of the 
nucleon[35-37, 48-53]; however, the difference in quark helicities could be attributed, 
in part, also to spontaneous chiral symmetry breaking[38]. In this connection, we 
think that the correct origin and the basic mechanisms for producing the Sivers 
asymmetry should be investigated more deeply. 

The initial and final state interactions may be described by the so-called link 
operator, introduced in the definition of bilocal functions in order to assure gauge 
invariance[25, 38, 39, 40, 41]. Moreover they cause also a nonvanishing [44, 45, 46, 
47]: in a scalar diquark model, this function is equal to [42] (see also ref.[47]). In 
the mechanism which generates quark transversity in an unpolarized nucleon, angular 
momentum conservation implies a change by one or two units of orbital angular 
momentum of the quark; this change can be connected to a pseudovector particle 
exchange, while the above mentioned initial or final state interactions are interpreted 
as Regge (or absorptive) cuts [44, 45, 46, 47]. 

From the above discussion it follows that quark-gluon interactions are essential 
for producing T-odd functions. Indeed, if such interactions are turned off, T-odd 
functions are forbidden by time reversal invariance[22] in transverse momentum space. 
On the contrary, they are allowed in the impact parameter space [48-53]: the Sivers 

7 



asymmetry can be viewed as a left-right asymmetry with respect to the nucleon spin 
in that space, where final state interactions produce a chromodynamic lensing for the 
struck quark[48-53]. 

As shown in appendix, T-odd functions can be related[97, 78] to the Qiu-Sterman 
[98, 99, 100] effect, which takes into account quark-gluon-quark correlations (see also 
ref.[101, 102]). This relation will be illustrated especially at the end of sect. 5, in 
connection with singly polarized Drell-Yan: T-odd functions turn out to produce an 
asymmetry whose Q 2 dependence coincides with the one obtained by assuming one 
of the above mentioned correlations[103, 104]. T-odd functions can be approximately 
factorized[38] - up to a sign, according as to whether the functions are involved in 
SIDIS or in Drell-Yan[38, 97] - if condition (5) is fulfilled [42]; otherwise one is faced 
with serious difficulties as regards universality of the effect [105]. 

2.4 Fragmentation Correlator 

The fragmentation correlator can be parameterized analogously to $, see subsect. 
2.2. We have, in the case of quark fragmentation into a pion, 

A = A e + A , (12) 

where, at leading twist, the T-even part is given by 

A e ~ X -k ■ n' + Dfi'_ (13) 

and the T-odd part reads 

Here k is the four-momentum of the quark, k± = k — (k ■ n'_)n' + — {k • n' + )n'__ and 
n'± are a pair of lightlike vectors, defined analogously to n±, but such that the space 
component of n'_ is along the pion momentum. /!q is an energy scale analogous 
to /io- Lastly D and H^- are fragmentation functions, D is the usual one, chiral 
even, while - the Collins function [22] - is chiral odd. It is important to notice 
that the latter function is interaction dependent, as well as the T-odd distribution 
functions: indeed, it has been shown [106] that this function would vanish in absence 
of interactions among partons. 
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3 SIDIS and Drell-Yan Asymmetries 



Now we deduce the expressions of the asymmetries involved in the two reactions 
considered, according to the formalism introduced in the previous section (see also, 
e. g., refs.[25, 29, 107, 108] for SIDIS and ref.[109] for Drell-Yan). As regards SIDIS, 
we treat the cases where the initial lepton is either unpolarized or longitudinally 
polarized, while the nucleon target is either unpolarized or tranversely polarized. On 
the other hand, concerning Drell-Yan, we consider the situations where at most one 
of the two initial hadrons (typically a proton) is transversely polarized. For our 
aims, the most relevant kinematic variables are two azimuthal angles, denoted as (j) 
and 4>s- In the case of Drell-Yan they are the azimuthal angles, respectively, of the 
momentum of the positive lepton and of the spin of the initial polarized hadron, in a 
frame at rest in the center of mass of the final lepton pair. Different frames, related to 
one another by rotations, have been defined: while the x-axis is usually taken along 
q_i_, the z-axis may be taken along the beam momentum - Gottfried- Jackson(G J) 
frame[110] -, antiparallel to the target momentum - U-channel (UC) frame [1, 2, 3] 
-, or along the bisector of the beam momentum and of the direction opposite to 
the target momentum - Collins-Soper (CS) frame[lll]. We shall discuss the frame 
dependence of the asymmetry parameters later on. As far as SIDIS is concerned, 
(j) and (fis are respectively the azimuthal angles - defined in the Breit frame where 
the proton momentum is opposite to the photon momentum - of the final hadron 
momentum and of the target spin vector with respect to the production plane. 



3.1 SIDIS Asymmetries 

The doubly polarized SIDIS cross section, with a longitudinally polarized lepton and 
a tranversely polarized nucleon, may be written as a sum of 4 terms, i. e., 

(da\ { da\ f da\ I ' da\ ( da\ . . 

Here we have singled out the unpolarized (UU), the singly polarized - either with a 
transversely polarized target, (UT), or with a longitudinally polarized beam, (LU) - 
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and the doubly polarized (LT) contributions. According to the formalism introduced 
in sect. 2, we get, at leading twist approximation, 

£e^ + [/> s 20], (16) 

a 

]T e 2 a [Slsin((f) + <f> s ) + S a 2 sm{<P - fa) 

a 

S$sin(3(j) - <j) S ) + S2-sin2(j)], (17) 
0, (18) 

+ £>So«(0 -&,)]. (19) 

a 

Here we have expressed the cross section in units a 2 xz 2 s/Q A , where s is the overall 
cm. energy squared. Moreover, omitting the flavor indices of the functions involved, 
we have 



U 


= A(y)T[w Uo J u D], 


(20) 


Ui 


= -C(y)^F[ WUl ,hi,Ht\, 


(21) 


Si 


= C(y)\S ± \^F[w Sl ,h 1T ,Ht], 


(22) 


s 2 


= A(y)\S ± \^F[ws 2 ,ft T ,D], 


(23) 


s 3 


= C(y)\S ± \4^F[™s 3 ,hi T ,Ht], 


(24) 


s. 


= \C(y)^F[w S4 ,hi L ,Ht], 


(25) 


D 1 


= \\ e -E(y)F[w Dl ,g 1L ,D], 


(26) 


D 2 


= X e l-E(y)^T[w D2 ,g 1T ,D}. 

Z fJ,Q 


(27) 



We have denoted by Xe and Sj_ respectively the helicity of the initial lepton and the 
transverse component of the nucleon spin vector, with |Sj_| = sincps and A = cos4>s- 
Moreover 

A(y) = l-y+l/2y 2 , C(y) = l-y, E(y) = y(2-y), (28) 
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where y ~ q /£ and q is the four- momentum of the virtual photon, such that 
\q 2 \ = Q 2 . Lastly, T is a functional [29], 

T\w, f,D] = J d 2 p ± w(p ± , <u_)f(p±)D[z, z(q ± - p ± )], (29) 

w, f and D being, respectively, a weight function, a distribution function and a 
fragmentation function. As to the weight functions, we have 

w Uo = w Dl = 1, (30) 

= ws 4 = 2u-p±u-k±-k±-p±, (31) 

w Sl = u k ± , w 52 = u! Da = u- p ± , (32) 

w Ss = 4(u • p ± ) 2 u • k ± - 2u • p ± k ± • p ± - u • k ± p^. (33) 

Here we have set u = qj_/q±, p± = p±/q± and kj_ = (qj_ — p±)/q±. Notice that the 
first two terms of the cross section (17) correspond respectively to the Collins and 
Sivers asymmetry [22, 20, 21]. 

3.2 Weighted Asymmetries in SIDIS 

The weighted asymmetries are defined as 

Aw = (34) 

Here brackets denote integration of the weighted cross section over q± and over the 
azimuthal angles defined above. W is the weight function, consisting of the Fourier 
component we want to pick up [see eqs. (16) to (19)], times (q±/ fio) Ha (q±/ Ho) nb , 
where n a and n b are respectively the powers with which p^ and kj_ appear in the 
functions w [see eqs. (30) to (33)]. For instance, the weight function corresponding 
to the Collins asymmetry is Ws 1 = (q±/ f/,o)sin((f> + (fis)- 

3.3 Drell-Yan Asymmetries 

In the case of singly polarized Drell-Yan with a transversely polarized proton, we have 
(see also ref.[109]) 

f<] ° ' = E e aK a + ^ a cos20], (35) 
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dV ^ 

, / UT a 



^ e 2 a [S?sin{<t> + <j> s ) + S 2 a sin(<j) - <j> s ) 



da 

I UT 

+ S' 3 a sin(3(f) - (f) S )}. (36) 



Here we have adopted the same approximation as before and have expressed the cross 
section in units a 2 /3Q 2 . Moreover 

U' Q = A'iy^iwu^hJ.l (37) 

U[ = CW-^JTwuiMrt], (38) 

Si = -C'(y)^F[w Sl ,h 1T ,htl (39) 

S' 2 = A'{y) q -±T[w S2 JtTJi\, (40) 

^3 = -C'(y)^ 7 ^[w S3 ,hi T ,hfl (41) 

A\y) = 1/2 -y + y 2 , C\y) = y{l-y) (42) 



and 



y = 1/2(1 + cos9), (43) 



^ being the polar angle of the positive lepton in one of the frames (GJ, UC, CS) defined 
at the beginning of this section. /i and /i' are energy scales relative to the two initial 
hadrons in the Drell-Yan process. The change of sign of the T-odd functions with 
respect to SIDIS has been taken into account in the coefficients S[, S' 2 and S' 3 , as 
already discussed at the end of subsect. 2.3. 



4 Determining fiQ 

Here we derive the appropriate value of the parameter /i for sufficiently large Q 2 . 
To this end we expand (see appendix) the correlator in powers of the coupling and, 
by exploiting the Politzer[112] theorem on equations of motion (see also ref. [113]), 
we compare the zero order term and the first order correction respectively with the 
T-even and with the T-odd parameterizations, eqs. (7) and (8). We shall show that 
the two procedures lead to consistent results. 
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Figure 1: Behavior of the asymmetry parameter v vs the dimensionless parameter 
p = q±/Q. Data are taken from refs.[l, 2]: circles correspond to y/s = 16.2 GeV, 
squares to y/s = 19.1 GeV and triangles to y/s = 23.2 GeV. The best fit is made with 
formula (61), A = 1.177. 

4.1 Spin Density Matrix 

In appendix we show that, in the case of a transversely polarized nucleon, one has [114, 
115], in the limit of g — > 0, 



Here m is the rest mass of the quark, such that p 2 = m 2 , and S q is (up to a sign) the 
quark Pauli-Lubanski vector, defined so as to coincide, in the quark rest frame [116], 
with the Pauli-Lubanski vector S of the nucleon in its rest frame. Now we compare 
the various Dirac components of the density matrix (44) with those of the T-even 
correlator (7), taking into account relation (A. 21) between S and S q , and 




(44) 




(45) 
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As a result we get the following relations for a free, on-shell quark[115]: 

\ ± hi T = (1 - ei )Aj./i 1T , (46) 

A±#rr = (i-e 2 )X±hiT- (47) 

Here 

A ± = - p± ■ S/V, (48) 

moreover e 1 ~ m/P and e 2 ~ m/2V are the correction terms due to the quark mass, 
which is small for light flavors. The terms of order O [(m 2 + p 2 ± ) / V 2 } have been 
neglected. 

In order to determine /j,q, we observe that the functions g±r, hyr and hj~ T , involved 
in formulae (46) and (47), are twist 2, therefore they may be interpreted as quark 
densities. For example, gyr is the helicity density of a quark in a tranversely polarized 
nucleon. Therefore it is natural to fix /x in such a way that gyr and h^ T are normalized 
like hiT- This implies, neglecting the quark mass, 

X ± = X ± , (49) 

and, according to eqs. (10) and (48), 

fi = V = -j=p ■ n_. (50) 

4.2 First Order Correction in the Coupling 

We show in appendix that the first order correction in g of the correlator - denoted 
as $i in the following - is T-odd and corresponds to a quark-gluon-quark correlation. 
This confirms that T-odd functions vanish if we neglect quark-gluon interactions 
inside the hadron. Moreover the result binds us to compare $i to the parameterization 
(8) of the T-odd correlator. As proven in appendix, the comparison yields 

fx oc V, (51) 

consistent with eq. (50). Result (51) is a consequence of the Politzer theorem, of 
four-momentum conservation and of kinematics of one-gluon exchange. But the T- 
odd functions hj~ and may be interpreted as quark densities (see subsect. 2.3), 
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provided they are properly normalized. Therefore we adopt for them the same nor- 
malization as for the T-even density functions - e. g., gyr - which involve the vector 
?7_i_. This leads again to result (50), analogously to the case of noninteracting partons. 
Therefore we conclude that, for sufficiently large Q 2 , the energy scale /io has to be 
identified with V. This is true both for some (T-even) functions of the parameteri- 
zation (7) - which survive also in absence of quark-gluon interactions - and for the 
T-odd correlator (8), which, on the contrary, depends on such interactions. Thus the 
conjecture proposed at the end of subsect. 2.2 is confirmed. Incidentally, eqs. (46) 
and (47) involve just T-even functions, therefore, according to the results deduced in 
appendix, these equations are valid up to terms of 0(g 2 ) and are expected to hold 
down to reasonably small Q 2 . 

4.3 Remarks 

At this point some important remarks are in order. 

a) Of course, we expect /io to be modified by nonperturbative interactions: for ex- 
ample, in the case of the already cited quark-diquark model[35-37, 44-47] (see subsect. 
2.3), the interference term scales with Q 2 , in agreement with the assumption /i = M 
[79, 25]. However, the virtuality of the exchanged gluon increases proportionally to 
Q 2 [117, 118], so that, at increasing Q 2 , the gluon "sees" the single partons rather 
than the diquark as a whole and eq. (50) appears more appropriate. To summarize, 
if we take into account the intrinsic transverse momentum of quarks, we are faced 
with the normalization scale /i , which, for large Q 2 , is equal to V, while for smaller 
Q 2 (such that nonperturbative interactions are not negligible) it is of the order of the 
hadron mass, in accord with a phase space restriction. 

b) Since, as already observed, the parameter /i is encoded in the four-vector rj± = 
P±/ 'fi>o, result (50) implies that, for sufficiently large Q 2 , the transverse momentum 
has to be treated as an effective higher twist. This agrees with the observation by Qiu 
and Sterman[98, 99, 100] that, owing to gauge invariance, transverse momentum has 
to be paired with a transversely polarized gluon, which, through quark-gluon-quark 
correlations, gives rise to a higher twist contribution. But, as already stressed, and 
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as shown in appendix, T-odd functions may be viewed as correlations of this kind. 
On the other hand, also the function g 1T , T-even and classified as twist 2, results to 
be suppressed for large V. All that casts some doubts on the correlation between the 
twist of an operator and the V dependence of the corresponding coefficient [75, 76], 
when transverse momentum dependent functions are involved. Indeed, in this case, 
although the Dirac operator commutes with the Hamiltonian of a free quark [95, 96], its 
mean value over the nucleon state may be suppressed, if the density function describes 
interference between two quark states with different orbital angular momenta. This 
occurs for those "soft" functions which are washed out by pj_-integration. 

c) It is worth comparing our approach to Kotzinian's[73], who starts from the 
approximate expression of the density matrix for a free ultra-relativistic fermion and 
adapts it to the case of a quark in the nucleon. He parameterizes the density matrix 
with the 6 twist-2, T-even functions that appear in the parameterization (7). Similar 
results are obtained by Ralston and Soper[79] and by Tangerman and Mulders [23]. 
The difference with our approach is that those authors do not take into account the 
Politzer theorem [112], which implies relations among the "soft" functions. 

4.4 Determining fiQ in the Fragmentation Correlator 

As regards the fragmentation correlator, we adapt our previous line of reasoning to 
the case of a quark fragmenting into a transversely polarized spin- 1/2 particle, say a 
A. For g — > one has 

A - p< = h$ + M')(D + H llb ?). (52) 

Here M' and S' are, respectively, the mass and the Pauli-Lubanski vector of the final 
hadron, whereas Hi is the transversely polarized fragmentation function; the other 
symbols are those introduced in subsect. 2.4. By comparing this limiting expression 
with a parameterization of A - analogous to eq. (7) as regards twist-2 terms - , we 
get (j% = k ■ n' + /V2. 
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Q (GeV) 



Figure 2: Behavior of the asymmetry parameter v vs the effective mass Q of the final 
lepton pair at fixed q±. y/s = 23.2 GeV Data from ref. [2], CS frame, and fitted with 
formula (61), A ■ q\ = 2.52 GeV 2 . 



5 Q 2 Dependence of Asymmetries 

Now we apply the result of the previous section to the processes considered in the 
present paper. It is convenient to take the space component of n_ along the direction 
of one of the two initial hadrons (Drell-Yan) or along the direction of the virtual 
photon (SIDIS). In both cases we get V ~ Q/2. Therefore we assume 

/*o = (53) 

the result being trivially extended to fi' Q and to /!q. As a consequence, we conclude 
that the azimuthal asymmetries illustrated in sect. 3 decrease with Q 2 . In particular, 
as regards SIDIS, we predict 

S 1 ,S 2 ,D 2 (xp, U 1 ccp 2 , S 3 <xp 3 (54) 

and 

n M c P" M 

K ~Q ' ~Q~' ( ' 
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Figure 3: Same as fig 2. = 19.1 GeV. A ■ q\ = 2.09 GeV 2 . 



where 



P = q±/Q- 



(56) 



Results (55) are consequences of the fact that A [see eq. (9)] is proportional to 
Q^ 1 for a transversely polarized nucleon. Such predictions might be checked by 
comparing data of experiments which have been realized (HERMES [4, 5, 6, 7] and 
COMPASS [9]) with those planned (CLAS[10]), which operate in different ranges of 
Q 2 . A strategy could be, for instance, to isolate the various Fourier components in 
the cross section [see eqs. (20) to (27)] by means of the weighted asymmetries and to 
study their Q 2 dependence. A particular remark is in order as regards the unpolarized 
SIDIS asymmetry, which we predict to decrease as l/Q 2 , just like the twist-4 cos20 
asymmetry arising as a consequence of the quark transverse momentum[119, 120, 73]. 
This makes the two asymmetries hardly distinguishable, but the last asymmetry can 
be parameterized, as well as the coscp asymmetry (the Cahn effect [34]), by means of 
the unpolarized quark density. 

Concerning Drell-Yan, the predictions are 



S[,S' 2 oc p, 





(57) 
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As regards U[, the result will be checked against unpolarized Drell-Yan data in the 
next section; the other three predictions could be verified, in principle, by comparison 
with data from experiments planned at various facilities, like RHIC[82], GSI[83, 84] 
and FNAL[85]. It is important to observe that, according to the approximation 
assumed in the present article, the asymmetry terms S[, S' 2 , S' 3 and U[ are invariant 
under rotations (see eqs. (38) to (41)) and therefore independent of the frame chosen 
(CS, GJ or UC). 



As a conclusion of this section it is worth recalling that the Drell-Yan single spin 
asymmetry, integrated over the transverse momentum of the final muon pair, was 
studied some years ago, in terms of a quark-gluon-quark correlation function, and it 
was found to decrease as Q' 1 [103, 104, 121, 122, 123] (see also refs. [98-102, 124]), 
consistent with our result (57). 
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Figure 5: Same as fig 2. Data from ref.[l], ^1 = 19.1 GeV. A ■ q\ = 2.71 GeV 2 . 

6 Azimuthal Asymmetry in Unpolarized Drell-Yan 

As is well-known, unpolarized Drell-Yan presents an azimuthal asymmetry. This has 
been seen, for example, in reactions of the type[l, 2, 3] 



(58) 



where N is an unpolarized tungsten or deuterium target, which scatters off a neg- 
ative pion beam. The Drell-Yan angular differential cross section is conventionally 

expressed as 



Ida 3 1 



, 1 + \cos 2 9 + nsin29cos(p H — vsin 2 9cos2(j)\ . 
a ak I 47r A + 3 V 2 / 



(59) 



Here ft — (0, <f>), and <fi being respectively the polar and azimuthal angle of the 
fi + momentum in the one of the frames defined in sect. 3. Moreover A, fi and v are 
parameters, which are functions of the overall center-of-mass energy squared s, of q 2 L) 
of Q and of the Feynman longitudinal fractional momentum xp of the muon pair with 
respect to the initial beam. 
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In the naive Drell-Yan model, where the parton transverse momentum and QCD 
corrections are neglected, one has A = 1, /i — v — 0. Therefore deviations of such 
parameters from the above predictions - observed experimentally both for A and u, 
while fi is consistent with [1, 2, 3] - can be attributed to transverse momentum or 
gluon effects, as illustrated in ref. [1]. The main contribution to the Drell-Yan cross 
section derives from QCD first order effects, typically from the partonic reactions 
qq — > #7* and qg — > qj* [1, 2, 3, 125], which also could account for the asymmetry 
parameter v [1, 2, 3, 126]. However, such effects fulfill the Lam-Tung[127] relation, 
which, instead, turns out to be rather strongly violated (see refs.[42, 1, 2, 3] and 
refs. therein). This fact has led people to propose alternative mechanisms [42, 43] 
for explaining the asymmetry. Furthermore we notice that data[l, 2, 3] exhibit for 
v a substantial independence of the frame chosen (CS, GJ, UC), just as predicted 
by T-odd functions (see the previous section), while first order perturbative QCD 
corrections would imply [126] a considerable frame dependence for that parameter. 

The behavior of Drell-Yan data may be understood by observing that the cross 
section is very sensitive also to power corrections [128, 129, 130] (see also refs. [131, 
132]). In particular, a A ^ 1 is obtained by assuming for reaction (58) a simple 
model of initial state interactions[129, 130], somewhat similar to the quark-gluon- 
quark correlations [98, 99, 100]. Here the Drell-Yan unpolarized cross section is of the 
type 

^oc|/o + A| 2 , (60) 

where f is the naive Drell-Yan amplitude and f\ consists of two terms (due to gauge 
invariance), describing one gluon exchange between the spectator quark of the meson 
and each active parton. It results [129, 130] |/o| 2 oc (1 — x) 2 (l + cos 2 8), \fi\ 2 oc p 2 cos 2 6 
and 29? /oA* oc (1 — x) psir^Ocos^cosipQ, where ip is the relative phase of the two 
amplitudes. This implies A — 1 oc p 2 /(l — x) 2 , in good agreement with datafl, 2, 3]. 
However, p depends crucially on ipo, moreover the third term of eq. (59) is absent. 
This could be recovered by inserting in eq. (60) a third amplitude, say /{, describing 
one gluon exchange between each active parton and the spectator partons of the 
nucleon: the missing asymmetry is reproduced by the interference term 2?R,f[f*, as 
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sketched at the end of subsect. 2.3, which amounts to recovering T-odd functions [42]. 

Indeed, the features of the parameters \i and v are suitably interpreted in terms 
of the correlator: comparison of eq. (35) with eq. (59) yields p, — and 



A ^ = A p 2 , (61) 



with 



A = ^< r h ^b} . (62) 

Here eqs. (37), (38), (42) and the second eq. (57) have been taken into account. We 
make some approximations concerning A . First of all, we neglect its q± dependence: 
for example, if we assume a gaussian behavior as regards the p_i_-dependence of the 
density functions, the q± dependence disappears in the ratio (62). Secondly we neglect 
the Q 2 evolution of the "soft" functions involved; such a dependence is expected to 
be quite smooth, as follows by assuming factorization and demanding factorization 
scale independence for the hadronic tensor [133]. Lastly, as told in subsect. 2.1, we 
neglect the Sudakov suppression: indeed, this effect, as well as the previous one, would 
imply a weak Q 2 dependence [92, 93], more complicated than the (approximate) p- 
dependence[92, 93] exhibited by data (see fig. 1 of the present paper and tables 1 to 
3 in ref.[l]). Therefore we approximate A by a constant. 

We fit formula (61) to the experimental results of v at different energies, both as 
a function of p (fig. 1) and as a function of Q at fixed q± (figs. 2 to 6), treating 
Aq as a free parameter. We stress that the p-dependence of v cannot be reproduced 
by the assumption p = M [23, 24, 25], not even taking into account the effects, 
just discussed, of QCD evolution and Sudakov suppression. This assumption would 
provide also a poor approximation to data of v versus Q at fixed q±. 



7 Conclusion 

We have studied the parameterization of the transverse momentum dependent quark 
correlator, both for distributions inside the hadron and for fragmentation processes. 
We are faced with the energy scale p , introduced in the parameterization for di- 
mensional reasons, and determining the normalization of some of the quark densities 
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(or fragmentation functions) involved. Comparison of the parameterization with the 
limiting expression of the correlator for noninteracting quarks yields fiQ = p-n_, con- 
trary to the usual [23, 24, 25] assumption, /j,q = M, which appears more appropriate 
for situations where nonperturbative interactions among partons are present. The two 
different assumptions lead to different predictions on the Q 2 dependence of azimuthal 
asymmetries in SIDIS and Drell-Yan. Our result agrees with previous approaches to 
azimuthal asymmetries, in particular with the Q 2 dependence predicted by quark- 
quark-gluon correlations[98, 99, 100] for Drell-Yan single spin asymmetry, and also 
with data of azimuthal asymmetry in unpolarized Drell-Yan. In particular, our in- 
terpretation of this asymmetry is considerably simpler than the one which could be 
obtained with the usual assumption about Hq. Further challenges to the two different 
theoretical predictions could come from future Drell-Yan experiments [82, 83, 84, 85] 
and from comparison between present [4, 5, 6, 7, 9] and incoming[10] SIDIS data. 
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Appendix 

We illustrate some features of the correlator; in particular, we expand this quantity 
in powers of the coupling and study in detail the zero and first order term of the 
expansion. The correlator is defined as [25] 



Here N is a normalization constant, to be determined below, and P, S) is defined 
in such a way that 



ip being the quark (or antiquark) field of a given flavor and \P,S) a state of a nucleon 
with a given four-momentum P and Pauli-Lubanski four- vector S, while p is the quark 
four-momentum. The color index has been omitted in ijj for the sake of simplicity. 
Moreover 



C{x) = Pexp [-igAj(x)] , with A T (x) = / XaA^dz", (A. 3) 



is the gauge link operator. Here g is the coupling and " P" denotes the path-ordered 
product along a given integration contour X, A a and A® being respectively the Gell- 
Mann matrices and the gluon fields. The link operator depends on the choice of 
X, which has to be fixed so as to make a physical sense. According to previous 
treatments [38, 25], we define two different contours, X±, as sets of three pieces of 
straight lines, from the origin to x loo = (±oo, 0, 0j_), from x loo to x^oo = x + , xj_) 

and from X200 to x = (x~ , x + , xj_); here the ± sign has to be chosen, according as 
to whether final or initial state interactions [38, 25] are involved in the reaction. We 
have adopted a frame - to be used throughout this appendix - whose z-axis is taken 
along the nucleon momentum, with x ± = l/\/2(£ ± z). 

T-even and T-odd correlator 




(A. 1) 




(A. 2) 




We set 



E(0) 



(A. 4) 
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where <3>± corresponds to the contour 2± in eqs. (A. 3), while &' E and <3>' select 
respectively the T-even and the T-odd "soft" functions. These two correlators contain 
respectively the link operators Ce{x) and £o(x), where 

Ce(0)(x) = ^P{exp -igA x+ (x) ± exp -igA x _(x)]} (A. 5) 

and Ax ± (x) are defined by the second eq. (A. 3). It is convenient to consider an 
axial gauge with antisymmetric boundary conditions [25], to be named G-gauge in 
the following. This yields 

At_(x) = -A i+ (x) (A. 6) 

and therefore 

Ce(x) = Pcos [gAj + (x) , jCo(x) = —iPsin ^gAj + (x) . (A. 7) 

Then the T-even (T-odd) part of the correlator consists of a series of even (odd) 
powers of g, each term being endowed with an even (odd) number of gluon legs. 
Moreover eq. (A. 6) implies that the T-even functions are independent of the contour 
(2+ or X_), while T-odd ones change sign according as to whether they are generated 
by initial or final state interactions [38]. In this sense, such functions are not strictly 
universal [38]. 

The two conclusions above, as well as the power expansion in the coupling, turn 
out to be gauge independent, since the the correlator is by definition gauge indepen- 
dent for any value of g and the same is true for any term in the expansion. As a 
consequence, the zero order term is T-even, while the first order correction is T-odd. 
This confirms that no T-odd terms occur without interactions among partons, as 
claimed also by other authors[35, 36, 37, 38]. 

Let us consider the expansion of $' in powers of g, i. e., 

$' = %-ig& 1 + ... 1 (A. 8) 



with 



[%h = I 7^ e ^< P > s\^(o)M^)\p, S) (A. 9) 
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and 

(< ^ = / ^^jr^d^(P,5|^.(0)A a A-(^)^(x)|P,5>. (A. 10) 

We stress that the term (A. 10) consists of a quark-gluon-quark correlation, analogous 
to the one introduced by Efremov and Teryaev[101, 102] and by Qiu and Sterman[98, 
99, 100]. Inserting expansion (A. 8) into eq. (A. 1), we get 

$ = $ + $i + (A. 11) 

where 

$ = N J & (p; P, S)dp~ and $i = -igN J $i(p; P, S)dp~. (A. 12) 

From now on we shall consider a transversely polarized nucleon. Then, according to 
our previous considerations, we may identify, at the two lowest orders in g, <3> with 
$ e and $i with $ , where $ e and $ G are given, respectively, by eqs. (7) and (8) in 
the text. Now we study in detail these two terms of the expansion. 



Zero order term 

We apply the Politzer theorem on equations of motion[112], i.e., 

(P, S$(0)£(x)(ip - m)ip(x)\P, S) = 0. (A. 13) 

Here = <9 M — igXaA 1 ^ is the covariant derivative. The result (A. 13) survives 
renormalization and applies also to off-shell quarks. Expanding C(x) in powers of g, 
at zero order the theorem implies that the quark can be treated as if it were on shell 
(see also ref.[113]). Then we consider the Fourier expansion of the unrenormalized 
field of an on-shell quark, i. e., 

* {x) = / wk^M 5 ( p ~ - ^f^) e "*"I><(p)'.(p)> ( A - 14 ) 

where m is the rest mass of the quark, s = ±1/2 its spin component along the nucleon 
polarization in the quark rest frame, u its four-spinor, c the destruction operator for 
the flavor considered and V is defined by the first eq. (10) in the text: in our frame 
V = p + 1 ' As regards the normalization of u s and c s , we assume 

«A = 2m, (P,5|4(p0cy(p)|P, 1 S) = (27r) 3 5 3 (p , -p)5 ss ,g s (p), (A. 15) 
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where p = (p + , p±) and q s (p) is the probability density to find a quark with spin 
component s along the nucleon (transverse) polarization and four-momentum p = 
(p~, p), with p" = \Jm 2 + p]_/2p + . For an antiquark the definition is quite analogous. 
Substituting eq. (A. 14) and the second eq. (A. 15) into eq. (A. 9) and into the first 
eq. (A. 12), we get 

(®oh = ^E?.(p)[«.(p)]i[ 5 «(p)L' ( A - 16 ) 

But [u s (p)]i[tis(p)]_; is nothing but the matrix element pij of the spin density matrix 
of the quark. Therefore, taking into account the first eq. (A. 15), eq. (A. 16) yields 

N 1 

$ o = ^E9.(P)^ + m )( 1 + 2 ^)- ( A - 1? ) 

where 2sS q is the Pauli-Lubanski vector of the quark in a transversely polarized 
nucleon. We normalize the correlator according to the definition (4) of the hadronic 
tensor, that is, demanding that it reduce to the spin density matrix in the limit of 
g ^ 0. Therefore 

N = 2V. (A. 18) 

Eq. (A. 17) can be conveniently rewritten as 

$ = i(^ + m)[/ 1 ( a :,pi)+75^iT(x,pi)], (A. 19) 

where we have set, according to the definitions of the density functions [in the scaling 
limit, q s (p) -> q s (x, p\)], 

fi= E Qs, h 1T = E 2sq s . (A. 20) 

s=±l/2 s=±l/2 

Eq. (A. 19) corresponds to formula (44) in the text. According to the Politzer 
theorem, renormalization modifies the functions /i and h±T, but not the structure of 
this expression. 

Now we express S q as a function of S. The two vectors do not coincide, since the 
spin operator for a massive particle has to be defined in the particle rest frame [116]. 
Taking into account the proper Lorentz boosts, we get, for a transversely polarized 
nucleon, 

S i = S + J ± ^ + 0(f ± ), (A. 21) 

m 

27 



with = p±/V and Aj_ = — S ■ r] ± . 



First order correction 

Now we consider the term (A. 10), which, as shown before, is T-odd. This term 
gives rise to a final (or initial) state interaction[103, 30, 97, 121, 35, 36, 37], in the 
sense that the spectator partons of a given hadron may exchange a gluon either with 
the final active quark (in SIDIS) or with the initial active quark of another hadron 
(in Drell-Yan). This kind of interaction selects the direction of the momentum of the 
gluon in the triple correlation (A. 10). Indeed, if the gluon is emitted by the spectator 
partons, it must have the same direction as the active quarks in that correlation; if 
absorbed, it must have opposite direction. This observation is quite important, as we 
shall see. 

We apply again the Politzer theorem, eq. (A. 13), now considering the first order 
correction of C(x). We get, adopting the G-gauge, 



(^-m)$;=.M, (A. 22) 

where 

Mij = I Wf^^ S ^ ^ - ^M*)\P, S) (A. 23) 

and ^ = \ a ^ a (x), ^2 = A a 4 a (^2oo)- The matrix M. is a quark-gluon-quark correla- 
tor, with the dimensions of a momentum; this matrix is a "soft" quantity, therefore 
independent of the "hard" scale p + . Eq. (A. 22) implies 

*; = 4±I>, (A. 24) 

with p 2 ^ m 2 . On the other hand, by considering the Fourier expansion of the gluon 
field in eq. (A. 23), and by applying again Politzer's theorem, we conclude that the 
other quark in the triple correlation is on shell. Now we show that 

p 2 oc G»+) 2 (A. 25) 

for p + — > oo, so that $' x decreases like {p + )^ 1 in that limit. 
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Our previous considerations and four-momentum conservation imply that the 
quarks in the triple correlation have four-momenta p and p ± k respectively, such 
that 

(p±k) 2 = m 2 . (A. 26) 

Here k is the four-momentum of the gluon in the triple correlation, the ± sign referring 
respectively to gluon emission and absorption. Since the gluon is exchanged between 
two color charges, it is space-like, moreover 

fc=(fc ,k), with < k ~ \k x \ ~ \k y \ « \k z \ = 0(p + ). (A. 27) 

But, according to the previous discussion, the kinematics of one gluon exchange 
demands k z to be positive for emission and negative for absorption. Then condition 
(A. 25) follows from eqs. (A. 27) and (A. 26). 

However, in order to justify the parameterization (8) for $ (p) in the text, one has 
to make an approximation. Indeed, the hadronic tensor in SIDIS and in Drell-Yan is 
not rigorously factorizable into two terms, if we include the mechanism of one gluon 
exchange illustrated before. As discussed by Collins in ref.[38] and papers therein, 
an approximate factorization may be assumed, provided the transverse momentum of 
the final hadron (in SIDIS) or of the final pair (in Drell-Yan) is much smaller than the 
"hard" scale, see condition (5). Under such a condition, taking into account results 
(A. 24) and (A. 25) and the T-odd character of this term may be parameterized 
in the following way: 

-ig<S>[ = ^{f&rP*^nW±Sl + ik\\t±,fo]}- (A- 28) 

Here K is a numerical constant and f£r and are two "soft" functions. Inserting 
eq. (A. 28) into the second eq. (A. 12) yields a parameterization for $ 1; which, by 
comparison to eq. (8), leads us to conclude that /io oc V. As discussed in subsect. 
4.2, a suitable choice of the normalization of the two T-odd "soft" functions (which 
uniquely fixes the constant K) leads to eq. (50), as in the noninteracting case. 
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